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We present a study of B -> X(3872)is: with X(3872) decaying to D*°D° using a sample of 657 
million BB pairs recorded at the T(45*) resonance with the Belle detector at the KEKB asymmetric- 
energy e + e~ collider. Both D*° — > D 0r y and D*° — > D°tv° decay modes are used. We find a peak 
of 50.lt";? events with a mass of (3872.9to;4 t°i) MeV/c 2 , a width of (3.9t?;|t?;?) MeV/c 2 and a 
product branching fraction B(B -> X(3872)K)xB(X(3872) -> D*°D°) = (0.80±0.20±0.10) x 10~ 4 , 
where the first errors are statistical and the second ones are systematic. The significance of the 
signal is 6.4a. The difference between the fitted mass and the D*°D° threshold is calculated to 
be (l.lto'4 to'jj ) MeV/c 2 . We also obtain an upper limit on the product of branching fractions 
B(B -> Y(3940)K) X S(Y(3940) -> D*°D°) of 0.67 x 10~ 4 at 90% CL. 

PACS numbers: 14.40.Rt, 13.25.Hw, 12.39.Mk 



The X(3872) was discovered by Belle in 2003 in 
B ± -> J/ipir+TT-^ [H with a mass of (3872.0 ± 
0.6 ± 0.5) MeV/c 2 , and was later confirmed by CDF 0, 
D0 and BaBar 0. It is one of the many new and 
unexpected hidden-charm states states recently discov- 
ered with masses around 4 GeV/c 2 . So far it remains 
unclassified; it docs not seem to be a pure cc char- 
monium state, but may be a D*D deuson d, a 
tetraquark state Q or a charmonium-gluon hybrid Q. 
The current average mass in the J/tjjTr + Tr~ channel is 
(3871.50 ±0.19) MeV/c 2 @. 

An important feature of the A (3872) is that its 
mass is very close to the D*°D° threshold ((3871.81 ± 
0.36) MeWc 2 0). The X(3872) was also observed 
by Belle [11| as a near-threshold enhancement in the 
D°D°ir a invariant mass spectrum of the B — > D a D°ir°K 
channel, with a peak at (3875.2 ± 0.7±°{% ± 0.8) MeV/c 2 , 
where the third error is from the uncertainty on the D° 
mass [10], a Gaussian width of (2.42 ± 0.55) MeV/c 2 , 
and a branching fraction B(B -> D°D°ir°K) = (1.22 ± 
0.3lt°;3o) x 10~ 4 . This initial study did not distinguish 
between decays via the D*°, and more general D°ir° fi- 
nal states. Looking in addition for the D*° — > Z?°7 decay 
is crucial to demonstrate the presence of A (3872) decay 
through a D*°. 

The BaBar collaboration recently published an ob- 
servation of the decay B -> A(3872)(-+ D*°D°)K 
with a 4.9c significance [HI]. The observed mass is 
(3875. lt^ ± 0.5) MeV/c 2 and the width is (3.0ti;| ± 
0.9) MeV/c 2 , with a product branching fraction B(B + — > 
X(3872)is:+) x B(X(3872) -> Z?* /) ) = (1.67 ± 0.36 ± 
0.47) x 10~ 4 . In the BaBar analysis, D*° candidates were 
subjected to a mass-constrained fit. 

Another new particle called A (3940) was discovered 



by Belle in the e+e~ -> DD*°,J/vb process [I3|. A 
state with the same mass, the F(3940) (also known 
as A (3945) [loll), was discovered by Belle in B — > 
ljJ/vjK [14|] and was later confirmed by BaBar, al- 
beit with a smaller mass [Hj]. The possibility that the 
A(3940) and F(3940) are the same state has not yet been 
ruled out. 

In this paper we measure the A (3872) — » D*°D° decay 
mode, followed either by D*° — > D°j or D*° — > Z?°7r°, in 
charged and neutral -B — > A(3872)A decays. Inclusion 
of charge-conjugate modes is implied throughout the pa- 
per. Furthermore, we use the notation D*°D° to indi- 
cate both D*"D° and D*°D°. The results are based on 
a 605 fb _1 data sample, corresponding to 657 x 10 6 BB 
pairs, collected at the Y(45) resonance with the Belle 
detector [li} at the KEKB asymmetric-energy e + e~ col- 



lider [17J , which includes the statistics used in the previ- 
ous Belle analysis 11 1. 

The Belle detector is a general purpose spectrometer 
with a 1.5 T magnetic field provided by a superconduct- 
ing solenoid. A silicon vertex detector and a 50-layer cen- 
tral drift chamber are used to measure the momenta of 
charged particles. Photons are detected in an electromag- 
netic calorimeter consisting of CsI(Tl) crystals. Particle 
identification likelihoods Ck and are derived from in- 
formation provided by an array of time-of-flight counters, 
an array of silica aerogel Cherenkov threshold counters 
and dE/dx measurements in the central drift chamber. 

Charged tracks are identified as kaons using a require- 
ment on the likelihood ratio Ck / (Ck + C n ), which has an 
efficiency of 88% for kaons and 10% for pions. Similarly, 
the charged pion selection has an efficiency of 98% for pi- 
ons and 12% for kaons. ir° candidates are reconstructed 
from pairs of photons with energies E 7 > 50 MeV and 
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with invariant mass in the range 118 MeV/c 2 < M 77 < 
150 MeV/c 2 . A mass-constrained fit is applied to ob- 
tain the four-momenta of the w° candidates. Kg candi- 
dates are reconstructed in the Kg — > tt + t:~ mode with 
the requirement \M„ n — m K o\ < 15 MeV/c 2 , where m K o 
is the Kg mass [13] ■ Requirements on the Kg vertex 
displacement from the interaction point and on the dif- 
ference between vertex and Kg flight directions are ap- 
plied. The Kg selection criteria are described in detail 
elsewhere 18]. A mass- and vertex-constrained fit is ap- 
plied to improve the four-momenta measurements of the 
Kg candidates. 
D° 



mesons are reconstructed in the modes D° — > 
K--K+, K-ir+ir Q , A^Tr+vr+vr-, Kgir+ir- and K~K+. 
The D° candidates are required to have a mass within 
±14 MeV/c 2 (±26 MeV/c 2 for the K-jr+n" mode) of 
the D° mass, 1864.8 MeV/c 2 [13]. This mass window 
width corresponds to ±4cr (±3cr for A~7r + 7r°). Mass- 
and vertex-constrained fits are applied to improve the res- 
olution of the four-momenta of D° candidates. D° candi- 
dates are combined with a photon or a ir° to obtain D*° 
candidates. The photon candidate is required to have an 
energy in excess of 100 MeV and shower shape variables 
that are consistent with an electromagnetic shower; the 
ratio of energy deposition in the central 3x3 and 5x5 
crystals of the cluster is required to be greater than 0.8. 
A mass window of ±27.5 MeV/c 2 for the L>*° ->■ D°j 
channel and ±6 MeV/c 2 for the D*° ->■ D°tt° channel is 
applied (±3<r), and a mass-constrained fit is used to im- 
prove the four-momenta of the D*° candidates; the mass 
is constrained to 2007.0 MeV/c 2 [13]. 

B mesons are reconstructed by combining a D*°D° 
candidate with invariant mass Md-d < 4.0 GeV/c 2 and 
a charged or neutral kaon candidate. To further reduce 
the background, at least one D°(D°) is required to decay 
to K~7r + (K + ir~). The beam-energy constrained mass 



Afbc = 



F 2 

beam 



where P, is the momen- 



tum of the ith daughter of the B candidate in the e 
center-of-mass (CM) system, is required to be larger than 
5.2 GeV/c 2 . The energy difference AE = E B - E hca-m , 
where Eb is the CM energy of the B candidate and 
-Ebcam is the CM beam energy, is restricted to the range 
|A_E| < 25 MeV/c 2 . Continuum e+e~ — > qq background 
events (q = it, d, s, c) are suppressed by requiring the ra- 
tio of the second and zeroth Fox- Wolfram moments [13] 
to be smaller than 0.3. 

The average B candidate multiplicity per event is 2.3 
for D*° -> D°7 and 2.7 for D*° -> D°7r°. We select the 
candidate with the smallest value of the quantity 



AM D o 
cm d0 



/ AM go 



( A(M D ,o - M D o) 

\ 0-(M D ,o-M D o) 



/ AE 

\CTAE 



, (1) 



where Ax is the deviation of the measured quantity x 
from its expected value and <j x the uncertainty in its mea- 
surement obtained using a Monte Carlo (MC) method, 
and the last term applies to the D*° — * D Q ir° channel 
only. 

MC samples are produced using the EVTGEN [20j 
generator. The X(3872) mass distribution is generated 
according to a relativistic Breit-Wigner function 



where 



BW(m) 
T(m) 



fimT(m) 



(to 2 — n 2 ) 2 + /i 2 T(m) 2 
[1 p(m) 



(2) 



To 



TO p(p) 



p(m) = -^—\/{m 2 - (m D o + m D ,o) 2 ) 
Zm 

Xy/(m 2 - (m D o - m D ,o) 2 ), 

and n and Tq are the nominal mass and width of the 
resonance, respectively, and p(rn) is the momentum of 
one of the daughters in the rest frame of its parent. The 
term ml 1 (to) in the numerator of Eq. |j2| behaves like 
a phase-space function, giving a smooth rise near the 
D*°D° threshold. 

The response of the Belle detector is simulated using 
a GEANT3-based program [2l|. Since the A (3872) is 
very close to the D*°D° threshold, its mass resolution 
varies rapidly with the D*°D° mass (M D *g). This reso- 
lution was studied in detail using MC simulations. It is 
modeled as a single Gaussian with a width given by the 
function a\JM.n*n — mo (where a and mo are free param- 
eters) shown in Fig. [1] At 3872 MeV/c 2 , the resolution 
is 0.13 MeV/c 2 . 



> 
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3980 4000 



FIG. 1: X(3872) mass resolution as a function of the X(3872) 
mass in the D°7 channel, obtained from MC with a X(3872) 
mass spectrum generated for a continuous range of masses 
from threshold to 4.0 GeV/c 2 . Crosses are Gaussian reso- 
lutions for various generated D*°D° masses; the curve is the 
result of a fit with a square root function. Very similar results 
are obtained for the _D°7r° channel. 
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A two-dimensional unbinned extended maximum like- 
lihood fit to Mbc and Md* o is performed. The M^ c 
distribution is described by a single Gaussian function 
for the signal and an ARGUS function for the back- 
ground; the Mjyjj distribution is described by a rela- 
tivistic Breit-Wigner function convoluted with the mass- 
dependent Gaussian resolution for the signal and a square 
root function for the background. In the D°j channel, 
the signal function also includes a broad higher mass 
distribution, corresponding to D*°(D°'y)D events in- 
correctly reconstructed as D°D*°('~fD ). The shape and 
fraction of this contribution is determined from MC. In 
the D°tt° channel, however, the reflection shape is too 
similar to the signal one to be distinguished. This is 
taken into account as a systematic uncertainty. Addi- 
tional components are the y(3940) signal, described by a 
relativistic Breit-Wigner function, and the non-resonant 
B — ► D*DK background, which peaks in M^ c but not 
in Md* d and is therefore described by a single Gaussian 
function in Mbc and a square root function in Md*d- 
The fitting procedure has been extensively tested using 
toy MC samples. 

First, each D*° decay channel sample is fitted sepa- 
rately. The yield, mass and width of the A (3872) peak 
are free parameters of the fit, as well as the F(3940) 
yield and the number of background and non-resonant 
B -> D*DK events. The F(3940) mass and its width 
are fixed to the values of Ref. [1J|. The results of the 
fits are presented in Table HI Since all the results are 
consistent we proceed with a combined fit. 

We subsequently perform a simultaneous fit to both 
D*° modes where the mass and width of the signal func- 
tion are constrained to have the same values in both 
modes, but the ratio of the yields in the D°7 and -D°7r 
channels is left free. Table Q] shows the fit results. We 
obtain a yield ratio of N D o 1 /N D o^o = 1. 16^° 371 which 
is consistent with the value of 0.92 we expect from 
MC with no non-resonant D°D°ir° contribution. We 
then fix the ratios of D°^/D°tt signal yields for the 
A (3872) and F(3940) from MC studies assuming that 
the D°tt° and D°-f come from a D*°. This fit gives 
50. lixii? events with a statistical significance of 7.9a (see 
Fig. [2] and Table [TJ) . We compute the significance from 
— 2 In(£o/£max)) where £0 and £ max are the likelihood 
values returned by the fit with the signal yield fixed at 
zero and left free, respectively. This quantity should be 
distributed as x 2 ( n doi = 3), as three parameters are free 
for the signal. The results of the simultaneous fits are 
consistent with the results of the individual fits. The 
distributions for M^ c and A (3872) mass closer to the 
A(3872) signal region are presented in Fig. [3] 

Next we fit the B+ -> A(3872)A+ and B Q -> 
A(3872)Ag modes separately, fixing the A (3872) mass 
and width to the values obtained with the simultaneous 
fit. Table U shows the results of these fits. Assuming the 
B° — > A (3872) A transition rate to be equal to twice the 



B° — ► X(3872)Kg rate, we obtain a ratio of branching 
fractions 



B(B" -> A(3872)A°) 
B(B+ -> A(3872)A+) 



= 1.26 ±0.65 ±0.06, 



which is consistent with unity. In this ratio, most of the 
systematic uncertainties cancel out. The remaining un- 
certainties are MC statistics, particle identification and 
reconstruction efficiencies of the K + and Kg, which com- 
bine in quadrature to a 5% uncertainty. 

The systematic uncertainties for the mass, width and 
signal yield are estimated by varying the fixed param- 
eters of the simultaneous fit: the resolution function 
is scaled by factors 0.9 and 1.2 obtained from the in- 
clusive D*° mass resolution study; the ratio between 
D°7 and D°tt° yields is varied according to the error 
in the relative branching ratio between D*° — > D°-f and 
D*° — > D°n° [10] and the uncertainty in their relative 
reconstruction efficiency; the F(3940) mass and width 
are fixed instead to the values from Ref. [l5j : the yield 
of F(3940) is fixed to zero and to the upper limit of 
40 events; the shape of the combinatorial background 
is changed from a square root to an inverted ARGUS 
function. Additional uncertainties due to possible bi- 
ases of the fit parameters are obtained using a large 
ensemble of toy MC experiments generated for differ- 
ent values of A (3872) mass and width. We also con- 
sider systematic uncertainties due to the contribution 
from D*°D° — > D°D*° mis-reconstruction by chang- 
ing the shape of the miss-reconstruction component to 
one obtained from MC generated with different values of 
A (3872) mass and width. Using MC we verify that the 
crossfeeds between D°D°j and D a D°ir° do not produce 
peaks in the signal region. Another systematic uncer- 
tainty on the mass is due to the D*°D° threshold mass, 
2m D o + A(m D , ~m D o) = (3871.80 ±0.41) MeV/c 2 [l(|. 
The systematic uncertainties for the branching fraction 
are estimated from the following sources: the number of 
BB events, D° branching fractions, track finding efficien- 
cies, K/ir identification efficiency, 7 or tt° detection effi- 
ciency, K s reconstruction efficiency, limited MC statis- 
tics and variation of the fixed parameters of the fit. Ta- 
ble [ill summarizes the systematic uncertainties. We ob- 
tain total systematic errors of 1q 54 MeV/c 2 for the mass, 
l°;o7 MeV/c 2 for the width and ±12% for the branching- 
fraction measurement. The significance of the signal in- 
cluding systematic uncertainties decreases to 6.4a. 

In order to determine the reconstruction efficiency for 
the branching fraction measurement, MC samples of 
B+ -+ A(3872)H D*°D°)K+ and B° -► X(3872)(^ 
D*°D°)K^ events are generated for D*° — > and for 
D*° -> D°ir°. The A (3872) is generated with a mass of 
3872.5 MeV/c 2 and a width of 4.0 MeV/c 2 . The total 
MC efficiencies multiplied by subdecay branching frac- 
tions are presented in Table HI 
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M(D*D) (MeV/c^) 

D°tv° (right). The result of 



M(D*D) (MeV/c z ) 

FIG. 2: Distributions of M D , D mass for M bc > 5.27 GeV/c 2 for D*° -» D° 7 (left) and for D*° 
the simultaneous fit is shown by the superimposed lines. The points with error bars are data, the dotted curve is the signal, 
the dashed curve is the background, the dash-dotted curve (barely visible) is the F(3940) component, and the solid curve is 
the total fitting function. 




5.29 
M bc (GeV/c 2 ) 



3893.8 



M(D*D) (MeV/c^) 



FIG. 3: Distributions of M bc for Md-d < 3.88 GeV/c 2 (left) and of M d *d for M bc > 5.27 GeV/c 2 (right); the top row is for 
D*° — > Z)°7 and the bottom row for D*° — > D n . The result of the simultaneous fit is shown by the superimposed lines. The 
points with error bars are data, the dotted curve is the signal, the dashed curve is the background, the dash-dotted curve is the 
sum of the background and the B — > D* DK component, the dot-dot-dashed curve is the contribution from D°-D° reflections, 
and the solid curve is the total fitting function. 
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TABLE I: Summary of results: the fitted mass, width, and yield of the X(3872) peak, and the total reconstruction efficiency, 
branching fraction and statistical significance for the various fits described in the text. 



Sample 


M x (MeV/c 2 ) 


T (MeV/c 2 ) 


Yield 


exg 




tr-* /-art A. \ 

B (10 ) 


a 


D*° -> D°7 (XK+ and XK°) 


3873.4 ± 1.0 


4.2±f;| 


26.2±?;° 


4.56 x 10" 


-4 


0.87 ±0.28 ±0.10 


4Aa 


D*° D°tt° (XK+ and XK°) 


3872.8 ± 0.7 


3.iit: 5 


22.0+& 7 


4.93 x 10" 


-4 


0.68 ±0.26 ±0.09 


6.8a 


All (free D°y/D°-K° ratio) 


3872.9l^j 




50.6+_\ 4 d 


9.49 x 10" 


-4 


0.81 ±0.20 ±0.10 


7.9a 


All (fixed D°j/D°tv ratio) 


3872.9±o.4 


0.3_ 1 4 


50.1±^.! 


9.49 x 10" 


-4 


0.80 ±0.20 ±0.10 


7.9a 


B+ -> XK+ 


3872.9 (fixed) 


3.9 (fixed) 




8.17 x 10" 


-4 


0.77 ±0.16 ±0.10 


7.6a 


B° XK° 


3872.9 (fixed) 


3.9 (fixed) 




1.32 x 10" 


-4 


0.97 ±0.46 ±0.13 


2.8a 



TABLE 11: Sources of systematic uncertainties for the fit- 
ted mass and width (MeV/c 2 ) and for the B -> A(3872)/sT 
branching fraction (%). 



Source 


Mass 


Width 


B 


Resolution function 


±0.04 


±0.09 


±0.3 


D° 1 /D %° yields ratio 


±0.01 


±0.00 


±0.2 


y(3940) parameters 


±0.01 


+0.00 
-0.32 


+0.0 
-4.0 


Background shape 


±0.00 


+0.00 
-0.14 


+0.0 
-2.2 


Fit bias 


+0.05 
-0.30 


+0.15 
-1.00 


+5.0 
-0.0 


D°-D° reflections 


±0.02 


±0.11 


±0.5 


D° and D*° masses 


±0.41 






Number of BB events 






±1.4 


D° branching fractions 






±2.4 


Tracking efficiency 






±5.0 


Particle identification 






±4.0 


7 or 7T° reconstruction 






±7.3 


K% reconstruction 






±4.5 


MC statistics (efficiency) 






±1.4 


Total (quadrature) 


+0.41 
-0.54 


+0.22 
-1.07 


±12 



The branching fraction, assumed to be equal for 
charged and neutral B modes, is 



B(B -> X(3872)K) x 2?(X(3872) -> D*°D°) 
= (0.80 ±0.20 ±0.10) x 10~ 4 , 



(3) 



with a statistical significance of 8.0<r. Using a Flatte-like 
parameterization [231 ] . with g = 0.3 and f p — 0.007, we 
obtain E f = (-14.9 ± 2.0) MeV/c 2 , consistent with ex- 
pected E f = -11 MeV/c 2 [23], and a signal of 65 ± 12 
events with a statistical significance of 8.8<r. The data 
statistics are not sufficient to distinguish between other 
fitting functions suggested in the literature 24 1 . 

The peak mass in this mode, Mx, like that in the 
BaBar analysis [l2| . should not be directly compared to 
the mass of the peak seen in J/ipir + n~ , or to the mass of 
the AT(3872) state itself. In this analysis, to improve sig- 
nal/background separation, D°7 and D°tt° combinations 
are each subjected to a mass window selection, and then 
a mass-constrained fit. This procedure returns masses 
above D*°D° threshold by construction, so the distribu- 
tion (and the mass and width of the peak) should be 
interpreted accordingly. Efforts in this direction have al- 
ready appeared in the literature 24, [2^], addressing a 



for the sum 
D*°D° and 



where B{X{3872) -> D*°D°) stands 
of branching fractions for AT (3872) — > 
X(3872) -> D*°D°. 

We obtain a mass of (3872.9 + o;4to'.5) MeV/c 2 and a this with the upper limit (O we get 
width of (3.9t 2 jt?; 2 ) MeV/c 2 ._ The difference in mass 
between the peak and the D*°D° threshold is 



preliminary version of the results in this paper. It should 
be noted that the mass measurement in our earlier pa- 
per ll|, while lacking the mass constraint, is also non- 
trivial to interpret, due to the role of the D* width, and 
interference, in the decay amplitudes 24, 25j |. 

For the y(3940) state, the simultaneous fit yields 7 ± 
21 ± 4 signal events and we set an upper limit of 

B(B -> r(3940)if)xB(y(3940) D*°D°) < 0.67xl0~ 4 

(5) 

at 90% CL. B y a veraging the branching fractions of 
Refs and [H, we obtain B(B Y(39A0)K) x 
S(y(3940) ujJ/ip) = (0.51 ± 0.11) x 10~ 4 ; combining 



B(y(3940)->o;.7/V>) 



> 



8(y(3940)^D*°n°) 

0.71 at 90% CL, to be compared with the 90% CL lim- 
its from Ref. 13;], /3(AT(3940) -> uJ/ijj) < 0.26 and 



SM = M x - m D ,o - m D0 = (l.l±g;| ±g-|) MeV/c 2 . (4) 

As a cross-check, different shapes are used as a sig- 
nal function to fit the data. Using a non-relativistic 
Breit-Wigner function truncated at the threshold, we 
obtain a mass of (3873.4 + °;^) MeV/c 2 , a width of 
(4.3± 2 ;4) MeV/c 2 and a signal yield of 39.6lg;i events 



B(X(3940)->^J/V>) 



< 



B(X(3M0) -> D*°D°) > 0.45, thus B(x(3940) ^ D . 0DO) 
0.58 with more than 90% CL, an incompatibility that 
suggests that the A" (3940) and the F(3940) are different 
states. 

In summary, we find a significant near-threshold en- 
hancement in the D*°D° invariant mass spectrum in 
B — > D*°D°K decays. The significance of this en- 
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hancement including systematic uncertainties is 6.4<r; 
significant signals are seen in both D*° — > D°-f and 
D°ir° modes. Using a relativistic Breit-Wigner, we ob- 
tain a mass of (3872.9±^t°-j) MeV/c 2 and a width of 
(3.9li;|t?;i) MeV/c 2 . The difference between the fitted 
mass and the D*°D° threshold is (l.ll^toi) MeV/c 2 . 
Note that a D*° mass constraint has been used in this 
analysis; the fitted mass of the AT(3872) peak is 2.3a 
lower than the value obtained by BaBar [jjj], where 
a similar constraint was used. For the Y"(3940) state, 
we set an upper limit on the B(B — * U(3940)i\) x 
2?(Y(3940) — > D*°D°) branching fraction which suggests 
that the AT(3940) and the U(3940) are different states. 
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